Complete mitochondrial genome of the frillneck lizard (Chlamydosaurus kingii, Reptilia; Agamidae), another squamate with two control regions.

Using PCR, the complete mitochondrial genome was sequenced in three frillneck lizards (Chlamydosaurus kingii). The mitochondria spanned over 16,761bp. As in other vertebrates, two rRNA genes, 22 tRNA genes and 13 protein coding genes were identified. However, similar to some other squamate reptiles, two control regions (CRI and CRII) were identified, spanning 801 and 812 bp, respectively. Our results were compared with another Australian member of the family Agamidae, the bearded dragon (Pogana vitticeps). The overall base composition of the light-strand sequence largely mirrored that observed in P vitticeps. Furthermore, similar to P. vitticeps, we observed an insertion 801 bp long between the ND5 and ND6 genes. However, in contrast to P vitticeps we did not observe a conserved sequence block III region. Based on a comparison among the three frillneck lizards, we also present data on the proportion of variable sites within the major mitochondrial regions.

Characterisation of the complete mitochondrial genome and 13 microsatellite loci through next-generation sequencing for the New Caledonian spider-ant Leptomyrmex pallens

The complete mitochondrial genome and a set of polymorphic microsatellite markers were identified by 454 pyrosequencing (1/16th of a plate) for the New Caledonian rainforest spider-ant Leptomyrmex pallens. De novo genome assembly recovered the entire mitochondrial genome with mean coverage of 8.9-fold (range 1-27). The mitogenome consists of 15,591 base pairs including 13 protein-coding genes, 2 ribosomal subunit genes, 22 transfer RNAs, and a non-coding AT-rich region. The genome arrangement is typical of insect taxa and very similar to the only other published ant mitogenome from the Solenopsis genus, with the main differences consisting of translocations and inversions of tRNAs. A total of 13 polymorphic loci were also characterized using 41 individuals from a single population in the Aoupinié region, corresponding to workers from 21 nests and 16 foraging workers. We observed moderate genetic variation across most loci (mean number of alleles per locus = 4.50; mean expected heterozygosity = 0.53) with evidence of only two loci deviating significantly from Hardy-Weinberg equilibrium due to null alleles. Marker independence was confirmed with tests for linkage disequilibrium. Most loci cross amplified for three additional Leptomyrmex species. The annotation of the mitogenome and characterization of microsatellite markers will provide useful tools for assessing the colony structure, population genetic patterns, and dispersal strategy of L. pallens in the context of rainforest fragmentation in New Caledonia. Furthermore, this paper confirms a recent line of evidence that comprehensive mitochondrial data can be obtained relatively easily from small next-generation sequencing analyses. Greater synthesis of next-generation sequencing data will play a significant role in expanding the taxonomic representation of mitochondrial genome sequences.

Microsatellite loci and the complete mitochondrial DNA sequence characterized through next generation sequencing and de novo genome assembly for the critically endangered orange-bellied parrot, Neophema chrysogaster

A suite of polymorphic microsatellite markers and the complete mitochondrial genome sequence was developed by next generation sequencing (NGS) for the critically endangered orange-bellied parrot, Neophema chrysogaster. A total of 14 polymorphic loci were identified and characterized using DNA extractions representing 40 individuals from Melaleuca, Tasmania, sampled in 2002. We observed moderate genetic variation across most loci (mean number of alleles per locus = 2.79; mean expected heterozygosity = 0.53) with no evidence of individual loci deviating significantly from Hardy-Weinberg equilibrium. Marker independence was confirmed with tests for linkage disequilibrium, and analyses indicated no evidence of null alleles across loci. De novo and reference-based genome assemblies performed using MIRA were used to assemble the N. chrysogaster mitochondrial genome sequence with mean coverage of 116-fold (range 89 to 142-fold). The mitochondrial genome consists of 18,034 base pairs, and a typical metazoan mitochondrial gene content consisting of 13 protein-coding genes, 2 ribosomal subunit genes, 22 transfer RNAs, and a single large non-coding region (control region). The arrangement of mitochondrial genes is also typical of Avian taxa. The annotation of the mitochondrial genome and the characterization of 14 microsatellite markers provide a valuable resource for future genetic monitoring of wild and captive N. chrysogaster populations. As found previously, NGS provides a rapid, low cost and reliable method for polymorphic nuclear genetic marker development and determining complete mitochondrial genome sequences when only a fraction of a genome is sequenced.

Amplification of multiple copies of mitochondrial Cytochrome b gene fragments in the Australian freshwater crayfish, Cherax destructor Clark (Parastacidae: Decapoda)

Non-coding copies of fragments of the mitochondrial genome translocated to the nucleus or pseudogenes are being found with increasing frequency in a diversity of organisms. As part of a study to evaluate the utility of a range of mitochondrial gene regions for population genetic and systematic studies of the Australian freshwater crayfish, Cherax destructor (the yabby), we report the first detection of Cytochrome b (Cyt b) pseudogenes in crustaceans. We amplified and sequenced fragments of the mitochondrial Cyt b gene from 14 individuals of C. destructor using polymerase chain reaction (PCR) with primers designed from conserved regions of Penaeus monodon and Drosophila melanogaster mitochondrial genomes. The phylogenetic tree produced from the amplified fragments using these primers showed a very different topology to the trees obtained from sequences from three other mitochondrial genes, suggesting one or more nuclear pseudogenes have been amplified. Supporting this conclusion, two highly divergent sequences were isolated from each of two single individuals, and a 2 base pair (bp) deletion in one sequence was observed. There was no evidence to support inadvertent amplification of parasite DNA or contamination of samples from other sources. These results add to other recent observations of pseudogenes suggesting the frequent transfer of mitochondrial DNA (mtDNA) genes to the nucleus and reinforces the necessity of great care in interpreting PCR-generated Cyt b sequences used in population or evolutionary studies in freshwater crayfish and crustaceans more generally.

Mitochondrial DNA sequence and gene organization in the Australian Blacklip Abalone Haliotis rubra (Leach)

The complete mitochondrial DNA of the blacklip abalone Haliotis rubra (Gastropoda: Mollusca) was cloned and 16,907 base pairs were sequenced. The sequence represents an estimated 99.85% of the mitochondrial genome, and contains 2 ribosomal RNA, 22 transfer RNA, and 13 protein-coding genes found in other metazoan mtDNA. An AT tandem repeat and a possible C-rich domain within the putative control region could not be fully sequenced. The H. rubra mtDNA gene order is novel for mollusks, separated from the black chiton Katharina tunicata by the individual translocations of 3 tRNAs. Compared with other mtDNA regions, sequences from the ATP8, NAD2, NAD4L, NAD6, and 12S rRNA genes, as well as the control region, are the most variable among representatives from Mollusca, Arthropoda, and Rhynchonelliformea, with similar mtDNA arrangements to H. rubra. These sequences are being evaluated as genetic markers within commercially important Haliotis species, and some applications and considerations for their use are discussed.

Complete mitochondrial DNA sequences of the decapod crustaceans pseudocarcinus gigas (Menippidae) and macrobrachium rosenbergii (Palaemonidae)

The complete mitochondrial DNA sequence was determined for the Australian giant crab Pseudocarcinns gigas (Crustacea: Decapoda: Menippidae) and the giant freshwater shrimp Macrobrachium rosenbergii (Crustacea: Decapoda: Palaemonidae). The Pse gigas and Mrosenbergii mitochondrial genomes are circular molecules, 15,515 and 15,772 bp in length, respectively, and have the same gene composition as found in other metazoans. The gene arrangement of M. rosenbergii corresponds with that of the presumed ancestral arthropod gene order, represented by Limulus polyphemus, except for the position of the tRNA^Leu(UUR) gene. The Pse. gigas gene arrangement corresponds exactly with that reported for another brachyuran, Portunus trituberculatus, and differs from the M. rosenbergii gene order by only the position of the tRNA^His gene. Given the relative positions of intergenic nonoding nucleotides, the “duplication/random loss” model appears to be the most plausible mechanism for the translocation of this gene. These data represent the first caridean and only the second brachyuran complete mtDNA sequences, and a source of information that will facilitate surveys of intraspecific variation within these commercially important decapod species.

Phylogenetic relationships among nine scallop species (Bivalvia : Pectinidae) inferred from nucleotide sequences of one mitochondrial and three nuclear gene regions

Current knowledge of the evolutionary relationships among scallop species (Mollusca: Bivalvia: Pectinidae) in the Indo-Pacific region is rather scanty. To enhance the understanding of the relationships within this group, phylogenies of nine species of scallops with the majority from coastal regions of Thailand, were reconstructed by maximum parsimony, maximum likelihood, and Bayesian methods using sequences of the 16S rRNA of the mitochondrial genome, and a fragment containing the ITS1, 5.8S and ITS2 genes of the nuclear DNA. The trees that resulted from the three methods of analysis were topologically identical, however, gained different levels of support at some nodes. Nine species were clustered into two major clades, corresponding to two subfamilies (Pectininae and Chlamydinae) of the three currently recognized subfamilies within Pectinidae. Overall, the relationships reported herein are mostly in accordance with the previous molecular studies that used sequences of the mtDNA cytochrome oxidase subunit I, and the classification system based on microsculpture of shell features and morphological characteristics of juveniles. Levels of divergences were different among genes (i.e., the 5.8S gene showed the lowest levels of nucleotide divergence at all levels, whereas the 16S rRNA showed the highest level of variation within species, and ITS2 gene revealed the highest level of divergence at higher levels).

Selection on mitochondrial variants occurs between and within individuals in an expanding invasion

Mitochondria are critical for life, yet their underlying evolutionary biology is poorly understood. In particular, little is known about interaction between two levels of evolution: between individuals and within individuals (competition between cells, mitochondria or mitochondrial DNA molecules). Rapid evolution is suspected to occur frequently in mitochondrial DNA, whose maternal inheritance predisposes advantageous mutations to sweep rapidly though populations. Rapid evolution is also predicted in response to changed selection regimes after species invasion or removal of pathogens or competitors. Here, using empirical and simulated data from a model invasive bird species, we provide the first demonstration of rapid selection on the mitochondrial genome within individuals in the wild. Further, we show differences in mitochondrial DNA copy number associated with competing genetic variants, which may provide a mechanism for selection. We provide evidence for three rarely documented phenomena: selection associated with mitochondrial DNA abundance, selection on the mitochondrial control region, and contemporary selection during invasion.

Diverse eukaryotes have retained mitochondrial homologues of the bacterial division protein FtsZ

Mitochondrial fission requires the division of both the inner and outer mitochondrial membranes. Dynamin-related proteins operate in division of the outer membrane of probably all mitochondria, and also that of chloroplasts – organelles that have a bacterial origin like mitochondria. How the inner mitochondrial membrane divides is less well established. Homologues of the major bacterial division protein, FtsZ, are known to reside inside mitochondria of the chromophyte alga Mallomonas, a red alga, and the slime mould Dictyostelium discoideum, where these proteins are likely to act in division of the organelle. Mitochondrial FtsZ is, however, absent from the genomes of higher eukaryotes (animals, fungi, and plants), even though FtsZs are known to be essential for the division of probably all chloroplasts. To begin to understand why higher eukaryotes have lost mitochondrial FtsZ, we have sampled various diverse protists to determine which groups have retained the gene. Database searches and degenerate PCR uncovered genes for likely mitochondrial FtsZs from the glaucocystophyte Cyanophora paradoxa, the oomycete Phytophthora infestans, two haptophyte algae, and two diatoms – one being Thalassiosira pseudonana, the draft genome of which is now available. From Thalassiosira we also identified two chloroplast FtsZs, one of which appears to be undergoing a C-terminal shortening that may be common to many organellar FtsZs. Our data indicate that many protists still employ the FtsZ-based ancestral mitochondrial division mechanism, and that mitochondrial FtsZ has been lost numerous times in the evolution of eukaryotes.

Complete mitochondrial DNA sequence of the Australian freshwater crayfish, Cherax destructor (Crustacea: Decapoda: Parastacidae): a novel gene order revealed

The complete mitochondrial DNA sequence was determined for the Australian freshwater crayfish Cherax destructor (Crustacea: Decapoda: Parastacidae). The 15,895-bp genome is circular with the same gene composition as that found in other metazoans. However, we report a novel gene arrangement with respect to the putative arthropod ancestral gene order and all other arthropod mitochondrial genomes sequenced to date. It is apparent that 11 genes have been translocated (ND1, ND4, ND4L, Cyt b, srRNA, and tRNAs Ser(UGA), Leu(CUN), Ile, Cys, Pro, and Val), two of which have also undergone inversions (tRNAs Pro and Val). The ‘duplication/random loss’ mechanism is a plausible model for the observed translocations, while ‘intramitochondrial recombination’ may account for the gene inversions. In addition, the arrangement of rRNA genes is incompatible with current mitochondrial transcription models, and suggests that a different transcription mechanism may operate in C. destructor.

Extracting DNA from museum bird eggs, and whole genome amplification of archive DNA

We present a comprehensive protocol for extracting DNA from egg membranes and other internal debris recovered from the interior of blown museum bird eggs. A variety of commercially available DNA extraction methods were found to be applicable. DNA sequencing of polymerase chain reaction (PCR) products for a 176-bp fragment of mitochondrial DNA was successful for most egg samples (> 78%) even though the amount of DNA extracted (mean = 14.71 ± 4.55 ng/µL) was significantly less than that obtained for bird skin samples (mean = 67.88 ± 4.77 ng/µL). For PCR and sequencing of snipe (Gallinago) DNA, we provide eight new primers for the ‘DNA barcode’ region of COI mtDNA. In various combinations, the primers target a range of PCR products sized from 72 bp to the full ‘barcode’ of 751 bp. Not all possible combinations were tested with archive snipe DNA, but we found a significantly better success rate of PCR amplification for a shorter 176-bp target compared with a larger 288-bp fragment (67% vs. 39%). Finally, we explored the feasibility of whole genome amplification (WGA) for extending the use of archive DNA in PCR and sequencing applications. Of two WGA approaches, a PCR-based method was found to be able to amplify whole genomic DNA from archive skins and eggs from museum bird collections. After WGA, significantly more archive egg samples produced visible PCR products on agarose (56.9% before WGA vs. 79.0% after WGA). However, overall sequencing success did not improve significantly (78.8% compared with 83.0%).